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Abstract 

The differential cross section for the dilepton production is calculated including Fermi motion of 
hadron constituents as well as emission from the ladders in the formalism of unintegrated parton 
distributions. We use unintegrated parton distributions which fulfil Kwiecihski evolution equa- 
tions. Both zeroth- and first-order (for matrix element) contributions are included. We calculate 
azimuthal angular correlations between charged leptons and deviations from the pt{l + ) = Pt(J~) 
relation. We concentrate on the distribution in dilepton-pair transverse momentum. We find inci- 
dent energy and virtuality dependence of the distribution in transverse momentum of the lepton 
pair. We study also azimuthal correlations between jet and dilepton pair and correlation in the 
(Pit(jet),P2t(l + l~)) space. The results are compared with experimental data of the R209 and UA1 
collaborations. 
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I. INTRODUCTION 



The Drell-Yan dilepton production is one of representative examples for which QCD 
collinear perturbative calculation can be performed order-by-order. Usually inclusive distri- 
butions are discussed. The most often studied observables are: da/dM ee or da/dxp, where 
M ee is invariant mass of the dilepton pair and Xp is the Feynman variable of the pair. In this 
paper we concentrate rather on one- and two-dimensional distributions which are singular 
in the standard collinear approximation. 

In the 0-th order collinear approximation (quark-antiquark annihilation) the transverse 
momentum of the dilepton pair (sum of transverse momenta of opposite sign charged lep- 
tons) is zero due to momentum conservation. Then the Oth-order result is not included in 
calculating the distribution in dilepton transverse momentum. The lowest nonzero contribu- 
tions are 1-st order quark-antiquark annihilation and QCD Compton. Typical for collinear 
approach they show singularity at small dilepton transverse momentum. 

Due to inter-quark interactions the quarks/ ant iquarks, constituents of hadrons, are not at 
rest and posses nonzero transverse momenta. Already this effect causes that the 0-th order 
process contributes to the finite transverse momenta of the lepton pair. Furthermore the 
emissions of gluons before the qq — > hard process causes an extra fc r smearing which, via 
momentum conservation, lead to finite transverse momenta of the dilepton pair (see Fig{T]) . 
The initial transverse momenta are often modelled effectively in terms of phenomenological 
Gaussian distributions [l|, |2| . The effect of Fermi motion as well as emission from the ladders 
can be easily included in the formalism of Kwiecinski unintegrated parton distributions 0] . 

In the present paper we wish to calculate differential cross section for dilepton production 
in the formalism of unintegrated parton distributions. We shall include both 0-th order and 
1-st order contributions. We shall concentrate on the distributions in dilepton transverse 
momentum. This observable is extremely sensitive in the 0-th order to the initial transverse 
momenta of partons. The transverse momentum of Drell-Yan pair was calculated within 
next-to-leading order perturbative QCD [H as well as in the resummation formalism in the 
impact parameter space [5|. Our approach differs in details from those approaches. 

Our results will be compared with experimental data for elementary proton-proton or 
proton-antiproton scattering. We leave analysis of proton-nucleus scattering for a separate 
publication. 



II. FORMALISM 

A. 0-th order Drell-Yan cross section 



The differential cross section for the 0-th order contribution can be written as: 



dy 1 dy 2 d 2 p lt d 2 p 2 t J vr ix \Qtx 2 (xiX 2 s) 2 

5 2 (Ki t + «2t - Pit - fix) [Fq f (xi, >4ti Pf) Fq f ( x 2, i4v Vf) \M (qq ^ e+ e~)\ 2 

+J 7 q f (x 1 , kI v h f ) F qj (x 2 , K 2 2t , up) \M(qq e+e-)\ 2 } , 



(2.1) 



where ^(xi, n 2 t ) and J 7 i(x 2 ,K 2 , t ) are unintegrated quark/ ant iquark distributions in hadron 
hi and h 2 , respectively. 
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FIG. 1: The diagram for the 0-th order Drell-Yan dilepton prodution with initial emissions from 
the ladders. 

The longitudinal momentum fractions are evaluated in terms of final lepton rapidities 
and transverse momenta: 



where mt = yp 2 + m 2 is a so-called transverse mass. 

The delta function in Eq.( 12.1l) can be eliminated as e.g. in Refs. 0,0,1- 
Formally, if the following replacements 

Fifa, k 2 u , /4) -> x lPi { Xl , fi 2 F )5{K 2 lt ) , 
J r j (x 2 , n 2 2t , n 2 F ) -> x 2 pj(x 2 , /J> 2 F )5(i4t) 

are done one recovers standard text book formulae. 
B. 1-st order Drell-Yan cross section 

In the first order in a s there are two types of diagrams: QCD Compton and quark- 
antiquark annihilation. A typical diagrams for corresponding subprocesses are shown in 




(2.2) 
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FIG. 2: The subprocess diagrams for the 1-st order Drell-Yan dilepton prodution with initial 
emissions from the ladders. 

For example the multi-differential cross section for the QCD Compton can be written in 
terms of unintegrated quark/ ant iquark and gluon distributions as: 



da(hih2 — > 7*^0 ^-j f d 2 K lU d 2 K 2 t ' 



dy\dy2d 2 pitd?p2t ~* J vr tt 167r 2 (xiX2s) 



5 (ku + K 2 t - Pit - P2t) [F g (xi, K lt , n F ) F q/ (x 2 , n F ) \M(gq -> 7*g)| 5 

+F qf (xi, K 2 lt , Hp) Fg(x2, i4 t , h f ) \ m ((19 -> 7*?) I 2 ] 



(2.4) 



Similarly the cross section for the first-order quark- ant iquark annihilation associated with 
gluon emission can be written as 



da(hih 2 — > 7*X) y-^ /" d 2 /^ Gp/t 2t I 



dyidy2d 2 pi t d 2 p 2 t ~* J vr tt 167r 2 (xi£ 2 s) 



5 (k u + K 2t - Pit ~ P2t) [^iffa, K w fl F ) J 7 qf (x 2 ,K 2t , fi F ) \M (qq -> 7* gf)! 1 

+J" 9/ (xi,^,/i|) Fq f (x2, K\ t: H F ) l M (<?<? -^7*#)| 2 ] 



(2.5) 



The delta functions in Eq. (12.41) and Eq. (12.51) can be eliminated as e.g. in Refs.[6|, 0, IsJ. 

The cross section for the emission of the dilepton pair can be expressed in terms of the 
cross section for the emission of time-like photon written above times a probability of the 
transition of the virtual photon into dilepton pair as: 



da(hifi2 — > l + l jX) a e 



dM 2 3ttM 2 



da(h x h 2 -> -f*jX) , (2.6) 



where Mu is the dilepton invariant mass. Please note that M 2 { = Q 2 , where Q 2 is virtuality 
of the time-like photon. 

C. Nonperturbative region of small Mu 

The formalism presented up to here applies in the perturbative region when the dilepton 
invariant mass Mu is not too small. How to calculate Drell-Yan production for small invariant 
masses (Mu < 1 GeV) is not completely clear and this issue was not discussed in the 
literature. In the region of very small photon virtualities the standard formulae presented 
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in the section above do not apply directly and some modifications are necessary (the same 
is true for soft-gluon resummation method). The reason is twofold. First of all the parton 
distribution may not exist for very small scales. Secondly, there is a singularity when photon 
virtuality approaches 0. 

In the present analysis we shall use the following extrapolation procedure: 
Firstly, to be able to use the perturbative Kwiecihski parton distributions the factorization 
scale for UPDFs is taken as: 

/4 = Ml + Q 2 S , (2.7) 

instead of \j? F = Mf t used usually in calculating Drell-Yan cross section for large dilepton 
invariant masses. 

Secondly, to avoid singularities inherent in the matrix element (s — > 0) we use a simple 
replacement proposed in Ref. 0, @] : 



s + Q 2 s , 

12 



t - t - Q 2 J2 , (2.8) 
u — > u — Q 2 /2 . 

The parameter Q 2 is to be adjusted to experimental data. The above procedure allows 
to avoid singularity when Mu — » 0. Please note that such a replacement does not change 
s + t + u. 

The procedure described above is very similar in spirit to that used in calculating the 
"deep inelastic" structure function F 2 for small photon virtualities [l(J 11 



For illustration in FigJS] we present distributions in transverse momentum of one of the 
leptons: electron or positron (left panel) and distributions in transverse momentum of the 
pair of leptons (right panel) for different values of the parameter Q 2 . In this calculation 
we have used unintegrated parton distributions to be discussed in the next subsection. The 
results depend somewhat on the value of the parameter. In principle, the shift parameter 
Q 2 could be adjusted to experimental data for extremely low dilepton invariant masses. It 
is not clear a priori if such a procedure can be successful. 



D. Unintegrated parton distributions 

Due to its simplicity the Gaussian smearing of initial transverse momenta is a good ref- 
erence for other approaches. It allows to study phenomenologically the role of transverse 
momenta in several high-energy processes. We define a simple unintegrated parton distri- 
butions: 

jf auss (x,K 2 ,n F ) = xpf l (x,ii 2 F ) ■ f Gauss {n 2 ) , (2.9) 

where p^ oll (x, n 2 F ) are standard collinear (integrated) parton distribution {i = g,q,q) and 
fcaussi 1 ^ 2 ) is a Gaussian two-dimensional function: 

The UPDFs defined by Eq. (12.91) and (12.101) is normalized such that: 

J Ff aUSS (x,K 2 ,H 2 F ) dK 2 = X P f\x^ 2 F ) . (2.11) 
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FIG. 3: Distributions in transverse momentum of the electron or positron (left panel) and in the 
transverse momentum of the pair (right panel) in the nonperturbative region of Q 2 G (0,1) GeV 2 at 
W = 200 GeV and for different values of the parameter Q 2 . In this calculation Q 2 = 0.25 (dotted), 
0.5 (dashed), 0.75 (dash-dotted), 1.0 (solid) GeV 2 (from top to bottom). Here midrapidity electrons 
and positrons 2/1,2/2 S (-1,1) were selected for illustration. 



Kwieciriski has shown that the evolution equations for unintegrated parton distributions 
takes a particularly simple form in the variable conjugated to the parton transverse mo- 
mentum. In the impact-parameter space the Kwieciriski equation takes the following simple 
form 



df NS (x,b,^) a s (fi 2 ) f 1 



dfi 2 


27T/1 2 „ 




- f N s(x,b 


df s (x,b, /i 2 ) 


a s (n 2 ) 


dfj, 2 


27T/1 2 „ 




+ P q9 (z)fc 


df G {x,b, n 2 ) 


a s (fi 2 ) 



dzP qq (z) 



B(.:-.r)./ ((l (^-/r) 



dzl Q(z — x) Jo((l — z)fib) 



Pqq(z)fs (-,b,H 



X 



- [ zP qq( Z ) + zP gq( Z )] fs(%, ^ 



d\i 2 2-Kjj, 2 



dz<e(z-x)J ((l- z)nb) 



Pgq( z )fs (-,b,H 



X 



- [zP gg (z) + zP qg (z)]f G (x,b^ 2 ) 



We have introduced here the short-hand notation 



fNS — fu fui fd fd i 

fs = fu + fu + fd+fd + fs + fs 



(2.12) 



(2.13) 
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The unintegrated parton distributions in the impact factor representation are related to the 
familiar collinear distributions as follows 



/ fc (x,6 = 0,/i 2 ) = |p fc (x,/i 2 ). (2.14) 

On the other hand, the transverse momentum dependent UPDFs are related to the integrated 
parton distributions as 

poo 

xp k (x,fi 2 )= dn 2 t T k {x, /t 2 ,/i 2 ) • (2.15) 
Jo 

The two possible representations are interrelated via Fourier-Bessel transform 

/oo 
db bJ (K t b)f k (x,b,iJ?) , 

,00 ; ; (2-16) 

f k (x,b,fi 2 ) = / dK t K t Jo{Ktb)Fk{x,K 2 ,n 2 ) . 



The index k above numerates either gluons (k=0), quarks (k> 0) or antiquarks (k< 0). 

While physically JF fc (x, k 2 , jj?) should be positive, there is no obvious reason for such a 
limitation for f k {x, b, fi 2 ). 



In the following we use leading-order parton distributions from Ref.[12j as the initial 
condition for QCD evolution. The set of integro-differential equations in b-space was solved 
by the method based on the discretisation made with the help of the Chebyshev polynomials 
(see (sj). Then the unintegrated parton distributions were put on a grid in x, b and /1 2 and 
the grid was used in practical applications for Chebyshev interpolation. 

For the calculation of inclusive and coincidence cross section for the heavy (time-like) 
photon production (see next section) the parton distributions in momentum space are more 
useful. These calculation requires a time-consuming multi-dimensional integration. An ex- 
plicit calculation of the Kwiecihski UPDFs via Fourier transform for needed in the main 
calculation values of (x±, k 2 J and (X2, n^t) * s n °t possible. Therefore it becomes a neccessity 
to prepare auxiliary grids of the momentum-representation UPDFs before the actual calcu- 
lation of the cross sections. These grids are then used via a two-dimensional interpolation 
in the spaces (xi, n\ t ) and (x 2 , k 2 J associated with each of the two incoming partons. 



III. A COMMENT ON 6-SPACE RESUMMATION 

There is an alternative approach to calculate transverse momentum distribution of the 
dilepton pair. In the Collins-Soper-Sterman formalism [13J (see alsoj^J), known also as b- 
space resummation, the cross section differential in dilepton invariant mass (Mu), rapidity 
of the pair (y) and transverse momentum of the pair (pt+) can be written as 



d(j ! f 1 

d 2 bexp(ip t b)W{b,Mu,x 1 ,x 2 ) = — I dbbJ {pt+b)W{b,Mu,x 1 ,x 2 ) 



dM 2 dydpi (2vr) 2 J " v "" / v ' ' ' ' 2tt 



(3.1) 
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The integrand function W(b, M 2 t ,pt+, xi, x 2 ) for the Drell-Yan pair production is: 

47r 2 a 2 



W(b,M l j,p t+ ,xi,x 2 



V a 



+ 



X2 

A-K 2 a 2 

9Mf s 



1 d& 
6 



/a/hi I £l, 



/b//i 2 




(3.2) 



y{M 2 h p t+ ,x 1 ,x 2 ) . 



The first term is naturally called resummation term and the function Y gives a correction 
which is negligible for small p t+ . The decomposition (13.21) is not free of ambiguities, like 
matching condition etc. Furthermore it is not easy to assure that in each corner of the phase 
space (y, Mu, p t+ ) the resummation part of the cross section is positive. For p t+ <C Mu 
the resummation term is much bigger than the correction term. The formula (13.11) together 
with (13. 2p has a singularity when Mu — > due to the photon propagator. In Eq. fl3.2j) the 
exponent in the Sudakov-like form factors reads 



S(b,Q 2 



dfi 2 

<b) ^ 



log 



A(a s (^)) + B(a s (^ 2 )) 



(3.3) 



The coefficients A and B can be expanded in the series of a, [13 . 

In this approach longitudinal momentum fractions are calculated as a rule as for the 
collinear Oth-order kinematics: 



Xi 



exp(+y)-p, 



x 2 = exp(-y) 



Mu 



(3.4) 



Please note that this is slightly different than in our case, where the transverse momenta of 
leptons explicitly enter into corresponding formulae. The soft-gluon resummation formula 
could be corrected for finite p t+ . 

The coefficient functions Cfi can be expanded in terms of a s 



C fa = 6 fa 5(l-Ci) + ... 
C fb = 6^6(1 -&) + ... 



(3.5) 



where the dots represent higher order terms in a s . Limiting to the resummation term in 
(13. 2p and keeping only first terms in the expansion (a bit academic approximation used here 
for simplicity) one gets 



2^,2 



W res (b,Mlp t , Xl ,x 2 ) 



9Mp 



E 



b 2 



Qf x ? 



b 2 



Qf X! 



Qf x ? 



b 2 



(3.6) 



S 



Obviously b -C 1/A in order that the pQCD is at work. The region of large b is of nonper- 
turbative nature. In order to supplement to the whole b space an extrapolation is needed. 
Usually this is done with the help of the following prescription: 



^"•=(TTWC^ <6 ' (37) 

Above bmax is a free, a bit arbitrary parameter [\lb 2 max > /Zq, where //g is a minimal possible 
factorization scale). To cut off contributions of large impact factor an extra form factor 
F NP (b,x%,x 2 ,...) multiplying W res (b, Mu,p t+ ,xi,x 2 ) is introduced. Above dots represent a 
set of free parameters. These parameters must be found by fitting to experimental data. 
Typical, not small, uncertainties are shown e.g. in Fig.l of Ref. [B| (see also flo|). 

Our approach here was formulated in the momentum space. Let us limit (as in the fa- 
space resummation) to distributions in y, p t + and M 2 t . Then the cross section can be written 
approximately (!) as 



da _ ^ af Y f d 2 k lt d 2 k 2t 
dyd 2 p t +dM 2 l ^-f sx\X 2 J n it 



(fq f (x 1 ,k 2 t ,ix 2 F )fg f (x 2 ,kl t ,fx 2 F ) + fg f (x 1 ,k' 2 t ,fj 2 F )f qf (x 2 ,kl t ,fx 2 F )) S 2 [k u + k 2t - p t+ 



(3.8) 

where af Y = fey-^jf-. Analogously as for the W and Z boson production (see Ref.[l~6l|) 
the spectrum in y, p t + and Mf t can be written also in terms of 6-space unintegrated parton 
distributions. The corresponding cross section then reads: 



da 47ra s 



dyd 2 p t+ dM 2 9M 2 s vr' 
f q /i(x u b, i4>)fq/ 2 (x 2 , 6, fi%) + fq/i(x u b, fM%) f q/2 (x 2 , 6, fx 2 F ) 



(3.9) 



In most of the evalutions presented in this paper up = M 2 { is used. 

We wish to note very similar structure of the 6-space resummation formula (13 .ip with 
(13. 6p and our formula (13.91) . Asumming a factorizable form for the nonperturbative form 
fctctor ^ * 

F» P (Q, 6, x u x 2 ) = F? P (Q, 6, Xl ) ■ F q NP (Q, b, x 2 ) (3.10) 

one could define (Vspace resummation unintegrated (anti) quark distributions in the b-space 
as: 

f* f GR (x 1} 6, Q 2 ) = F q (b, Xl , Q 2 ) [x iqf (x u /i 2 (6)) + ...] exp (~S(b, Q 2 

} , ( (3.1D 

f q s f GR (x 2 ,b,Q 2 ) = F,(b,x 2 ,Q 2 ) [x 2 q f (x 2 ,fi 2 (b)) + ...] exp I --S(b,Q 2 



1 This is justified if the main nonperturbative effects are due to the merging of initial (anti)quarks in bound 
(nonpcrtubative) hadrons which causes that initial partons are not at rest and posses internal momenta. 
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FIG. 4: Distributions in azimuthal angle between electron and positron in proton-proton scattering 
for RHIC energy y/s = 200 GeV and different windows in Q 2 specified in the figure. The dashed 
lines are for UPDFs with Gaussian smearing (do = &o = 1 GeV -1 ) and the solid lines are for 

Kwiecihski UPDFs with b = 1 GeV -1 . Here \? F = Q 2 was chosen. No cuts on lepton transverse 
momenta were applied. 



The 6-space resummation is useful for calculating transverse momentum distribution of 
the dilepton pair, a variable Fourier-conjugated to b. For other correlation observables 
between e + and e~ exact kinematics and new kinematical variables must be used. This 
seems not possible in the framework of the 6-space resummation. Our approach, formulated 
in the momentum space and with explicit kinematics of each lepton, is better suited for 
calculating observables like or -5 — ^ discussed in the present paper. 



IV. RESULTS 



A. 0-th order component with transverse momenta 



In the standard collinear approach in 0-th order approximation leptons of opposite charge 
are produced back-to-back, i.e. the relative azimuthal angle between them is fixed and equals 
to 180°. The situation changes if one includes "initial" transverse momenta of partons which 
annihilate producing dilepton pair. This is shown in FigJU 

We have selected RHIC energy (yfs = 200 GeV) for illustration. Since the plot is purely 
theoretical we did not make any cuts on electron/positron transverse momenta and rapidi- 
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FIG. 5: The influence of cuts on lepton transverse momenta on distributions in azimuthal angle 
between electron and positron for = 38.8 GeV (E772 experiment). Here \j? f = Q 2 was chosen. 

ties. We show results for narrow bins in photon virtuality (= square of the dilepton invariant 
mass) specified in the figure. Of course, the cross section strongly depends on photon virtu- 
ality; the smaller virtuality the bigger cross section. This is mainly a kinematical effect. This 
can be seen by comparison of solid (Kwiecihski UPDFs) and dashed (Gaussian smearing) 
lines. In this calculation ctq of the Gaussian distribution was adjusted to bo in the Kwiecihski 
UPDFs. For small virtualities the two distributions almost coincide. At large virtuality they 
differ much more which can be understood as effect of evolution of UPDFs from the initial 
scale /1q to Q 2 . 

How the azimuthal angle correlations depend on cuts in lepton transverse momenta ? In 
Figj5]we show some examples for the E772 experiment at Fermilab. The cuts not only lower 
the cross section but also modify the shape of azimuthal correlations. 

In collinear Oth-order approximation the electron and positron transverse momenta com- 
pansate each other, i.e. pu{e + ) = Vn{ e ~)- m the ^-factorization approach this condition 
is relaxed. As an example in Figj6] we show two dimensional maps in the region of rela- 
tively small transverse momenta {pu,P2t < 5 GeV). In this calculation we did not put any 
constraint on photon virtuality (dimuon invariant mass). In the left panel we show results 
obtained with Eq. fl2.ll) . As discussed in the theoretical section real perturbative calcula- 
tion requires presence of large scales. Therefore to be precise the perturbative calculation 
at small transverse momenta is not reliable. This can be better understood by looking at 
FigJT] which shows correlations of transverse momenta and photon virtualities. In the right 
panel of Figj6]we show in addition calculation with shifted scales and kinematical variables. 
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FIG. 6: Distribution in (pit(e + ),P2t ( e ~)) f° r zeroth-order Drell-Yan in proton-proton collisions at 
yfs = 200 GeV. In the left panel standard procedure is used with /j, F = Q 2 and in the right panel 
in addition Q 2 S = 1 GeV 2 is used as described in subsection IIC. In this calculation -1 < 2/1,2/2 < 1. 




FIG. 7: Average value of the photon virtuality on the (pit(e + ),p2t(e )) plane for zeroth-order 
Drell-Yan in proton-proton collisions at yfs = 200 GeV. Here /fp = Q 2 and -1 < 2/1, 2/2 < 1- 

A careful inspection of both panels shows that the differences are only at small transverse 
momenta. A more orthodox approach would be to completely exclude the region of small 
transverse momenta of leptons, which means also small dilepton invariant masses. 

Very interesting observable, which is singular in collinear approximation in leading-order, 
is the distribution in the transverse momentum of the dilepton pair (jh+). In FigJS]we show 
such a dependence on the incident center-of-mass energy for two different bins in photon 
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FIG. 8: Distribution in pt+ for different beam energies specified in the figure and for two different 
windows in Q 2 . Here fi 2 F = Q 2 . 

virtuality. In general, the bigger energy the broder the distribution in p t+ 2 . The bigger 
energies correspond to smaller values of quark/ ant iquark longitudinal momentum fractions. 
The effect of broadening of the distribution is larger for larger photon-virtuality which is 
taken as a factorization scale in the Kwiecihski UPDFs. 

The effect of broadening is summarized in Fig 19] where is show average value of the lepton- 
pair transverse momentum. The difference between different bins in photon virtuality is 
due to QCD evolution effect (we use photon virtuality Q 2 as a factorization scale in the 
Kwiecihski unintegrated parton distributions). Similar effects were already discussed for 
dijet |7| and photon-jet correlations [8]. The two curves in Figj9] would coincide if the 
evolution effect would be neglected, as was done e.g. in Ref. 0]. 

We show the effect of broadening also for the RHIC energy y/s = 200 GeV where we 
show distribution in p t+ for different bins in photon virtuality. The effect of broadening 
is inherently related to the Kwiecihski UPDFs where the initial kf t and/or A;J* smearing 
depends on the factorization scale and on x\ and/or x 2 (see discussion in [1, 0, @|)- 

Finally we wish to confront our calculation with existing data for the Drell-Yan dilepton 
production. In Fig{10] we show our results with different values of the parameter 6 111 
the Kwiecihski UPDFs [3]. The parameter 6 quantifies nonperturbative effects which are 
beyond evolution effects embodied in the Kwiecihski evolution equations. The dominant 
effect is probably the Fermi motion of nucleon constituents. This effect is often neglected 
in the standard QCD calculations. We get quite good description of the R209 collaboration 
data already in the zeroth-order with bo — 1 - 2 GeV^ 1 . This is a bit surprising, as the zeroth- 
order contribution is usually neglected in orthodox collinear approach. It seems therefore 



2 In collinear approach the distribution would be delta function in transverse momentum of the pair. 
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FIG. 9: Average value of pt+ as a function of center-of-mass energy for two different windows of 
Q 2 . Here [i\ = Q 2 . 

indespensible to include higher-orders in the fc r factorization approach. We shall discuss this 
issue in the next section. 



B. 1-st order component with transverse momenta 

Let us now discuss contribution of processes of one order higher than in the previous 
section, with hard subprocesses shown in Figj2l These diagrams have to be inserted between 
two partonic ladders. 

The correlation in azimuthal angle between jet and the dilepton pair is shown in FigJTTl 
One can see a strong deviation from the collinear back-to-back kinematics caused by the 
transverse momenta of initial partons. The most top (thick) lines are for the full phase 
space. The intermediate lines are with extra cuts on e + e~ and jet rapidities. The most 
bottom lines are for extra cuts on transverse momenta of the associated jet Pt(jet) > 5 GeV. 
The latter cut changes drastically the shape of distributions which are peaked now more at 
= 180°. 

In Fig{T2] we show corresponding transverse momentum distribution of the dilepton pair 
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FIG. 10: Distribution in pt + for zero-order Drell-Yan in proton-proton collisions for W = 62 GeV. 
Different curves correspond to different values of the 6o parameter in the Kwiecihski UPDFs. The 
experimental data of the R209 collaboration are taken from 17]. 



for RHIC energy = 200 GeV. For comparison we show also zeroth-order contribution 
discussed in the previous section. The zeroth-order contribution dominates at small trans- 
verse momenta, while the first-order contribution at transverse momenta larger than about 
5 GeV. Here we include both Compton and annihilation processes. The first-order contribu- 
tions have no singularity at pt(e + e~) as in the collinear approach. The zeroth and first-order 
contributions overlap only in a very limited range of pt(e + e~) ~ 4 GeV. The zeroth-order 
fc r factorization component includes some higher order contributions via UPDFs. However, 
because the two contributions occupy rather different parts of the phase space no severe 
double counting is expected. 

When integrated the first-order contribution is significantly smaller than the zeroth-order 
one. This can be better seen in Fig JT3] where we show rapidity distributions. We show the 
zeroth-order (LO) and two /^-factorization Compton components: QCD Compton (qg — > 
qe + e~ and gq — > qe + e~) and quark-antiquark annihilation (qq — > ge + e~ and qq — > ge + e~). 

Let us present now the first-order contributions together with existing data. In Fig {14] we 
present separately two first-order contributions: Compton and quark-antiquark annihilation. 
The sum of the both contributions is shown by the thin solid line. This contribution is 
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FIG. 11: Distribution in azimuthal angle between jet and dilepton pair in proton-proton collisions 
for the RHIC energy y/s = 200 GeV for the QCD Compton (solid lines) and quark-antiquark 
annihilation (dashed lines). The results for -1 < y(e + e~),y(jet) < 1 are shown with thin lines and 
the results for -4 < y(e + e~), y(jet) < 4 are shown with thick lines. 



about factor of 4 smaller than the R209 collaboration data. For comparison we show also 
the zeroth-order contribution. The Oth-order contribution is much larger than the lst-order 
one. The situation may change at larger transverse momenta. The sum of the 0th- and lst- 
order well describe the transverse momentum data. The scenario discussed here seems quite 
different than the one presented in the text book by Field 19( where the data were discribed 
by a convolution of the collinear first-order component and some extra, somewhat arbitrary, 
Gaussian smearing and where the zeroth-order contribution was completely ignored. 

In Fig {15] we compare distributions obtained in collinear and fcj-factorization approach. 
The difference can be seen in small transverse momentum region. Above p t+ > 4 GeV 
the two approaches practically give the same results. The singularities seen for collinear 
approximation disappear when transverse momenta of initial partons are included. 

The first-order contribution becomes dominant at larger transverse momenta of the dilep- 
ton pairs. This is shown in Fig{16] where we present transverse momentum distribution of 
dimuons of opposite charge for proton-antiproton scattering at the center-of-mass energy 
y/s = 630 GeV. We confront results of our first-order fc r factorization calculation against 
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FIG. 12: Distribution in pj + for the RHIC energy y/s = 200 GeV. We show separately zero-order 
(dashed line) and first-order Compton (dotted line) contributions. Here -1 < y{e + e~) < 1 and 
rapidity of the jet for the first-order contribution is in the full phase space. 



UAl collaboration experimental data |18J which were measured at p t ^ > 10 GeV. We show 



our results for two different choices of the factorization scale: (a) [i 2 F = ML (left panel), 
(b) fj,jp = +p|„„ (right panel). The experimental data are well described within theo- 
retical uncertainties. In this case the zeroth-order contribution (not shown in the figure) is 
concentrated at pt lW1 < 10 GeV. 

In the first-order collinear calculations the transverse momentum of the dilepton pair is 
completely balanced by the transverse momentum of the associated jet (quark or antiquark 
for the Compton process and gluon for the quark-antiquark annihilation). This strict balance 
is not longer true when transverse momenta of initial partons are taken into account. In fact 
the disbalance can be a measure of the transverse momenta of the intial partons. In FigJTTl 
we show two-dimensional distributions in (pit(jet),p 2 t(e + e~)) for Compton (upper panels) 
and annihilation (lower panels) contributions for \/s = 200 GeV and a narrow window in the 
photon virtuality. We see broad distributions of the strength along diagonal pu = P2t with 
a smearing of the order of a few GeV. This smearing is a consequence of the convolution 
of two unintegrated parton distributions embodied in Eq. (12.4ft and ( 12.51) . The broadening 
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FIG. 13: Distribution in y(e + e~) for the RHIC energy -^/s = 200 GeV. We show separately zero- 
order (solid line) and first-order QCD Compton (dotted line) and quark-antiquark annihilation 
(dashed line) contributions. The rapidity of the associated jet for the first-order contributions is 
integrated in the full phase space. 

strongly depends on the choice of the factorization scale (compare left and right panels). 

In Figfl8]we show similar distributions for the UAl collaboration experiment. Here the 
range of transverse momenta of the jet and dilepton pair is much larger. For such a broad 
range of transverse momenta it looks as if the cross section is concentrated on the diagonal 
Pit{j e t) — P2*(A i+ A i ~)- A closer inspection shows a smearing similar as shown previously for 
proton-proton scattering at yfs = 200 GeV (see Fig JTTI) . 



V. CONCLUSIONS 

We have calculated both zeroth- and first-order contributions to dilepton production 
in the formalism with transverse momenta of initial partons taken into account. In these 
calculations we have used Kwiecihski unintegrated parton distributions which include both 
smearing in parton momentum due to nonperturbative effects in hadrons before collision as 
well as extra smearing due to QCD evolution effects in the collision process as encoded in 
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FIG. 14: Distribution in transverse momentum of the dilepton pair in proton-proton collisions. We 
show both first order contributions: QCD Compton (dotted line) and quark-antiquark annihilation 
(dashed line). For comparison we show also the zeroth order contribution (dash-dotted line). 

the Kwiecihski evolution equations. 

We have calculated correlations in azimuthal angle between both charged leptons as 
well as correlations in the two-dimensional space of transverse momentum of the positron 
and transverse momentum of the electron. Both effect of the Fermi motion and effect of 
subsequent emissions from the ladder lead to deviations from the delta function in rela- 
tive azimuthal angle centered at = 7r (collinear case) and deviations from p t (electron) = 
p t (positron) condition. The shape of the distribution in transverse momentum of the pair 
depends both on incident energy and virtuality of the time-like photon. This is a strightfor- 
ward consequence of the QCD evolution encoded in the Kwiecihski equations. We predict 
larger smearing in transverse momentum of the dilepton pair for larger dilepton masses. 
The existing experimental data at y/s = 62 GeV can be well explained by the zero-order 
component by adjusting the parameter responsible for nonperturbative effects of internal 
motion of partons in hadrons. This is rather in odds with the explanations in the literature, 
where the data are explained by an extra convolution of the first-order contribution with a 
Gaussian smearing function. In orthodox collinear the zero-order contribution is completely 
ignored. The smeared zeroth-order contribution discussed here may be partly responsible 
for missing strength in the spectrum of so-called nonphotonic electrons [20]. 

We have also calculated dilepton transverse momentum distribution in the first order for 
the matrix element. Inclusion of initial transverse momenta removes singularity at Ptj+i- = 
0. The first-order contribution dominates only at larger transverse momenta of the pair and 
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FIG. 15: Distribution in transverse momentum of the dilepton pair for the R209 collaboration 
experiment. The collinear (thin lines) distributions are compared to ^-factorization (thick lines) 
distributions for the QCD Compton (dotted) and quark-antiquark annihilation (dashed). 

is smaller than the zeroth-order contribution at low transverse momenta. The inclusion of 
initial transverse momenta leads naturally to decorellation of relative azimuthal angle of a 
jet and dilepton pair (in the first-order collinear approximation they are emitted back-to- 
back). We have also discussed analogous decorrelations on the (j>t{jet),p t (l + l~)) plane. The 
initial transverse momenta lead to sizeable deviations from the collinear condition pt(jet) = 
p t (l+l-). 

Finally we wish to make a comment on possible double counting. In principle, our leading 
order contribution contains diagrams which look like first order diagrams. The standard 
pQCD first-order result contains large log(Q 2 / pf) . For fixed Q 2 it happens when p\ is small. 
Fortunately it happens numerically that then (in our approach) the first-order result is much 
lower than the almost purely nonperturbative origin zeroth-order result. In principle, the 
double counting happens when the zeroth- and first-order results are comparable. As shown 
in our calculation this happens in a very narrow interval of p t+ . So we expect rather small 
double counting. 
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FIG. 16: Dilepton invariant mass distribution for first-order Drell-Yan processes in proton- 
antiproton scattering for two different scales: fi'p = Q 2 (left panel) and \xp = Q 2 + p 2 ^^ (right 
panel) at ^ = 630 GeV. Here Q 2 G (1,2. 5 2 ) GeV 2 and y 2 G (-1.7,1.7). The contributions of 
Compton (dotted) and annihilation (dashed) are shown separately. The experimental data of UA1 
collaboration are taken from 18]. 
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FIG. 17: Two-dimensional distributions in pu(jet) and P2t(& + & ) f° r the first-order Compton 
contributions (upper panels) and for the first-order quark-antiquark annihilation contributions 
(lower panels) in proton-proton collisions at yfs = 200 GeV and Q 2 £ (4.5,5.5) GeV 2 . Left panels 
are for = Q 2 and right panels for \j? F = Q 2 +p 2 ee . No cuts on rapidity of the jet as well on the 
rapidity of the dilepton pair were applied. 
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FIG. 18: (pit(jet),p2t (wO) distribution for first-order Drell-Yan processes in proton-antiproton 
collisions at ^/s = 630 GeV. Left panel for QCD Compton and the right panel for annihilation. 
Here [i\ = Q 2 and Q 2 = (1,2.5 2 ) GeV 2 and -1.7 < y{fi + >-) < 1.7. 
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